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Mitochondria are semi-autonomous organelles regulated by a complex network of proteins that are vital for 
many cellular functions. Because mitochondrial modulators can impact many aspects of cellular 
homeostasis, their identification and validation has proven challenging. It requires the measurement of 
multiple parameters in parallel to understand the exact nature of the changes induced by such compounds. 
We developed a platform of assays scoring for mitochondrial function in two complementary models 
systems, mammalian cells and C elegans. We first optimized cell culture conditions and established the 
mitochondrial signature of 1,200 FDA-approved drugs in liver cells. Using cell-based and C elegans assays, 
we further defined the metabolic effects of two pharmacological classes that emerged from our hit list, i.e. 
imidazoles and statins. We found that these two drug classes affect respiration through different and 
cholesterol-independent mechanisms in both models. Our screening strategy enabled us to unequivocally 
identify compounds that have toxic or beneficial effects on mitochondrial activity. Furthermore, the 
cross-species approach provided novel mechanistic insight and allowed early validation of hits that act on 
mitochondrial function. 



Mitochondrial dysfunction not only typifies a range of rare inherited mitochondrial diseases^'^, but it has 
become clear that defects in mitochondrial function are also associated with age-related diseasel The 
concept of maintaining or repairing mitochondrial function may hence offer new preventive and 
therapeutic strategies for several of these rare and common diseases that still require efficient treatments. 

Amongst the models that are available to explore mitochondrial defects and therapeutics, the mouse is the most 
widely used because of the existence of a powerful genetic toolkit and the physiological similarities to humans'*'^. 
For example, mitochondrial defects of COX- deficient mice can be reversed using the AMPK agonist AICAR^. 
Other compounds targeting similar pathways, i.e. resveratrol and nicotinamide riboside, restore mitochondrial 
function in mice fed with a high-fat diet^"^°. Notably, the beneficial effects of resveratrol also applied to human 
fatty acid oxidation- deficient fibroblasts^^ and obese patients^^. C. elegans is also used to study mitochondria, 
especially in the context of longevity^^"^^. The worm lives only weeks versus years for mice, while the ease of 
genetic modifications render the exploration of mitochondrial pathways accessible. Unfortunately, however, 
mouse and worm models are limited in throughput, precluding screening for new mitochondrial drugs^. One 
of the alternative approaches that has been developed in recent years involves phenotypic screening in both 
primary and immortalized cell lines^'^^"^^. These studies have shown that it is possible to identify compounds that 
impact on mitochondrial function without targeting a specific molecular pathway. Nevertheless, this type of 
approach requires a solid set of validated follow-up assays to elucidate how mitochondrial modulators work, and 
to narrow down several dozens of hits to compounds that are worth testing in animals. 

Here, we report the development of an assay platform composed of mouse liver cells (Hepal.6) and C. elegans 
for the systematic identification and characterization of compounds that improve mitochondrial metabolism 
(Figure 1). These assays are intended to assess compounds for their efficiency or toxicity on mitochondrial 
function. To validate our method, we first optimized cell culture conditions to ensure that the cells relied on 



SCIENTIFIC REPORTS j 4 : 5285 | DO!: 1 0.1 038/srep05285 



1 



B 



O 

E 



CO 

E 



c: 

CT3 
CD 



Redox potential 
ATP level 



Lactate production 

I 

Oxygen consumption 
(basal, uncoupled, acute) 



Effect on development 



Oxygen consumption 
(basal, uncoupled) 



In depth characterization 



Caloric restriction 
Exercise 
Resveratrol 

i 




Figure 1 | Assays platform for the identification and validation of novel mitochondrial modulators. (A) Workflow to assess mitochondrial function 
going from hundreds of compounds to the best candidates (top to bottom) in mammaUan cells and C. elegans. After hits have been identified during a 
high-throughput screening, they are assessed for their efficiency and toxicity on mitochondrial fijnction using this assay platform. Once the best 
candidates have been sorted out, they can be further investigated in both models, and ultimately in vivo. (B) The different pathways that are assessed and 
the respective markers that were measured (black boxes). This scheme is non- exhaustive and represents only the key regulatory steps that are surveyed. 
More details can be found in the main text. 



oxidative metabolism. We next identified compounds with a strong 
impact on mitochondrial function out of the —1,200 FDA-approved 
drugs in liver cells using a phenotypic strategy simflar to that 
developed previously^^. Two pharmacological classes of compounds 
emerged in our hit list, imidazoles and statins. We further character- 
ized the effect of these two pharmacological classes on mitochondrial 
function using our set of assays. In summary, our study reports the 
potential of the use of an intact organism, such as C. elegans, in 
combination with mammalian cells to identify and validate hits that 
affect mitochondrial function at a very early stage of the drug dis- 
covery process. 

Results 

Designing a multi-level screening platform in mammalian cells and 
C. elegans. We developed a platform of assays using mammalian cells 
and C. elegans to identify and validate compounds that modulate 
mitochondrial function (Figure lA-B). As a starting point, we 
adapted a phenotypic high-throughput screening strategy that was 
developed and initially validated in muscle cefls^^. During the 
screening and validation of mitochondrial modulators, the main 
objective is to rapidly establish the toxicity and efficacy of the 
compound. Several cell culture assays were employed for this 
purpose. For example, measurement of ATP levels is a robust way 
to assess toxicity^^'^° (Figure lA-B). Metabolic redox balance, such as 
NAD^/NADH ratio (Figure lA-B), can be measured indirectly with 
the use of fluorescent redox probes^\ Nicotinamide adenine dinucleo- 
tide (NAD^) is not only a cofactor involved in the redox reactions, but 
it is also a co -substrate of the sirtuin family of deacetylases, which 
control mitochondrial metabolism at various levels^^ (Figure IB). 
Fluorescent probes that accumulate inside mitochondria propor- 
tionally to mitochondrial membrane potential (A\|/M) and mitochon- 
drial abundance (e.g. tetramethylrhodamine or TMRM) can inform 
both about efficacy and toxicity of compounds^^ (Figure lA-B). Once 
hits are identified in this initial screen, the optimal dose for each 
compound can be determined in cells by using the same assays as 
used during high-throughput screening (Figure lA). Mitochondrial 
bioenergetics, or energy flow through mitochondria, is then assessed 
by monitoring oxygen consumption and lactate production^^'^^ 



(Figure lA-B). All these assays are transposable to almost every cell 
line, although they are representative only of the effect of a compound 
on the tissue from which the cell line is derived. 

In this regard, C. elegans offers the possibflity to score the impact of 
compounds on mitochondria and physiological homeostasis at the 
scale of an entire organism. A first and easy parameter to assess 
toxicity of a compound is to look at its impact on worm development 
(Figure lA). Indeed, mitochondrial mutants were shown to have 
either impaired or delayed larvae growth^^. To further characterize 
mitochondrial bioenergetics, oxygen consumption can also be mea- 
sured in worms^^ (Figure lA-B). 

In our model, we propose that only compounds with effects on 
mitochondria in both models are further validated in more expensive 
and time-consuming assays in both species (Figure lA). 

Optimization of the cell culture conditions. We developed our 
mammaUan assays in the mouse hepatocyte cell line Hepal.6, a cell 
line that is robust and easy to culture, and is derived from a tissue 
where mitochondria are abundant. We first optimized cell culture 
conditions to ensure that the cefls solicit mitochondria to produce 
energy. We monitored cellular energy utilization by comparing 
different nutrient sources, 1) high glucose (5 g/1, HG); or 2) low 
glucose (1 g/1) supplemented with 30 |iM oleic acid (LGO) 
(Figure 2A). While HG can be metabolized through glycolysis and 
mitochondrial oxidation, LGO forces the cells to actively rely on 
mitochondrial function. Indeed, LGO medium induced a marked 
decrease in two assays indicative of glycolysis, including redox poten- 
tial, which is sensitive to NAD^/NADH ratio, and extra- cellular 
acidification rate (ECAR), a marker of lactate production (Figure 2A), 
In contrast, ATP level and oxygen consumption rate (OCR) were both 
increased in LGO medium (Figure 2A). We also observed an increase in 
the NAD^/NADH ratio in Hepal.6 cells grown in LGO, in agreement 
with reduced redox potential (Figure 2B). These results show that LGO 
medium induces a metabolic shift from glycolysis to oxidative 
phosphorylation (OXPHOS). 

We next tested whether cefl confluence impacts on Hepal.6 meta- 
bolism by measuring OCR and ECAR at different levels of confiuence 
(Figure 2C). Whfle OCR increased with cefl number, ECAR 
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Figure 2 | Optimization of the cell culture conditions. ( A-D) Reduction in glucose concentration and supplementation with oleic acid shifts cells toward 
oxidative metabolism. (A) Redox potential, ATP level, oxygen consumption rate (OCR) and extra-cellular acidification rate (ECAR) measured in 
Hepal.6 cells after 24 hours incubation in the indicated growth medium. (B) NAD^/NADH ratio determined in Hepal.6 cells after 24 hours incubation 
with high glucose (5 g/l,HG); or low glucose (1 g/1) supplemented with 30 |aM oleic acid (LGO) (C) Hepal.6 cells shift toward oxidative metabolism with 
increasing cell confluence. OCR values are represented as a percentage of the value measured with 15,000 cells per well. ECAR values are represented as a 
percentage of the value measured with 500 cells per well. (D) Hepal.6 cells are only sensitive to oUgomycin at higher cell density (12,000 cells/well versus 
1,000 cells/well). ATP content was normalized for each cell density over the vehicle treated control cells. For A and B, Statistical significance was 
determined by Student's t-test, with * P<0.05, P<0.01, *** P<0.001. Graphs represent mean ± SEM. 



decreased until a certain threshold, suggesting that Hepal.6 cells shift 
from glycolysis to OXPHOS to produce ATP as confluence increases 
(Figure 2C). To confirm this observation, we treated Hepal.6 cells at 
different cell confluences for 24 hours with the mitochondrial ATP 
synthase inhibitor oligomycin, and measured residual ATP levels 
(Figure 2D). In line with our results on the metabolic switch at higher 
confluence, these cells were insensitive to high doses of oligomycin at 
low confluence, while becoming sensitive at higher confluence 
(Figure 2D). Based on these results, we decided to perform all assays 
when the cells are at 95% confluence and using LGO as cell culture 
medium. 

Screening of the Prestwick library of FDA-approved drugs in 
Hepal.6 cells. For the full validation of our high-throughput 
screen (HTS) we tested viability (ATP level), redox potential and 
mitochondrial membrane potential (A\|/m) in Hepal.6 cells 
exposed to a collection of —1,200 mostly FDA-approved drugs in 
the Prestwick library (Figure 3). We obtained good separation 
capacity for ATP and redox potential, with a Z-factor of 0.58 and 
0.85, respectively, while the A\|/m assay had a Z-factor of 0.01, 
suggesting that this assay could potentially lead to more false 
positives and negatives (Figure 3A). The absence of significant 
correlation between the different parameters shows that the assays 
are complementary and non-redundant (Figure 3B-C). 

Compounds were considered as hits when inducing a signal 
change superior to DMSO ± 2aDMso and when deviating from 
normal distribution, i.e. when squared distance between observed 
quantile and theoretical quantile was above a certain threshold 
(Figure 3D-E, methods section). 53 hits were identified after apply- 
ing this two-steps selection method (Figure 4A). Hierarchical clus- 
tering was used to identify the different profiles of the hits, i.e. the 
combinations of the scores in each of the assays (Figure 4A). We 
focused our subsequent validation work on pharmacological classes 
that were enriched in our screen by profiles. This strategy increased 
our chances to characterize an effect that is due to the common 
mechanism of action of a pharmacological class, rather than an 
off-target effect of an isolated compound. One cluster contained 
six hits that were particularly toxic to the cells, as all three parameters 
decreased (Figure 4A-B). It included the antineoplastic drug camp- 
tothecine (S, + ), and niclosamide, which has been reported to inhibit 
OXPHOS directly^^ (Figure 4A-B). A larger cluster encompassed 
compounds that increased redox potential, without affecting ATP 
or A\|/m (Figure 4 A and C). Amongst these compounds, five belong 
to the class of antifungal imidazoles (Figure 4 A and C). Imidazoles 
block demethylation of lanosterol through the inhibition of sterol 
14a-demethylase (CYP51A1)^^'^^, an essential step in the synthesis of 



ergosterol in fungi and cholesterol in mammals. Other reports show 
that their antifungal effect is also mediated by an induction of 
uncoupling and swelling of rat liver mitochondria^^. Imidazoles 
henceforth represent a class of compounds that are potentially toxic 
to the mitochondria. A third remarkable cluster regrouped com- 
pounds that induce A\|/m, without affecting ATP or redox potential 
(Figure 4A and D). Strikingly, this set of A\|/m inducers contained 
three statins (Figure 4 A and D). Statins inhibit one of the initial steps 
in cholesterol biosynthesis, but this also affects the synthesis of coen- 
zyme QIO (CoQlO), heme- A and substrates for protein prenyla- 
tion^^. Therefore, statins represent a class of compounds that may 
also induce mitochondrial stress, but have an overall demonstrated 
beneficial impact on human health and longevity^^. 

Considering the overrepresentation of statins and imidazoles in 
our hit list, we decided to further characterize their effect on mito- 
chondrial function using our mammalian and C. elegans assays 
(Figure lA). 

Impact of imidazoles on metabolism in cultured cells and C. 
elegans. We first performed dose-response tests for all our mam- 
malian assays using the four imidazoles identified, i.e. econazole, 
miconazole, tioconazole and bifonazole (Figure 5A-H). Tiocona- 
zole and bifonazole showed toxic effects on ATP levels, but only at 
30 jiM (Figure 5A). Miconazole and bifonazole increased redox 
potential by up to 50% at 10 and 30 |iM, while econazole and 
tioconazole gave milder responses (Figure 5B). All four imidazoles 
decreased A\|/m by at least 50% at 10 and 30 |iM (Figure 5C). The 
reason why these compounds did not decrease A\|/m significantly 
during the HTS (Figure 4A and C) is explained by the high variability 
of this assay, leading to a higher threshold of significance when 
testing multiple samples. All these changes were independent on 
cell number, as assessed by DNA quantification (Figure 5D). Long- 
term treatment of Hepal.6 (24 hours) cells with imidazoles induced 
a metabolic shift toward glycolysis, as ECAR was maintained 
(Figure 5E) while basal OCR significantly decreased (Figure 5F). 
Uncoupled OCR was also significantly reduced after imidazoles 
treatment, indicating a decrease in maximal respiratory capacity 
(Figure 5G). Measurement of OCR after short-term treatment 
showed that econazole, miconazole and bifonazole inhibit cellular 
respiration within a few minutes (Figure 5H), suggesting that 
imidazoles inhibit the mitochondria directly and acutely. This 
toxic effect of imidazoles was confirmed in C. elegans, as they all 
strongly disturbed worm development at 10 jiM and higher doses 
(Figure 51). 

The fact that C. elegans are cholesterol auxotrophs and are supple- 
mented with cholesterol during the experiments^^'^"'^^ underscores 
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Figure 3 | Screening of 1,200 FDA approved drugs in the Hepal.6 cell line. (A-E) Screening of the Prestwick hbrary of FDA approved drugs, using the 
three basal assays and oxidative cell culture conditions, i.e. low glucose concentration ( 1 g/1), supplementation with oleic acid at 30 [iM and cell density of 
95% the day of the assays. (A) Cell-based assays screening separation bandwidth. Plots represent the cumulative density of points over Z scores of the raw 
values for all the points of the screening, including negative (blue points) and positive controls (red points) and samples. The redox potential assay has the 
best separation capacity with narrow peaks of negative and positive controls on the density plot and a Z factor close to 1. (B-C) Correlations between the 
three cell-based assays show poor level of relationship between the different parameters. (B) Pearson's p correlation coefficients and their corresponding 
regression r^ value were calculated for normalized ratios of the three parameters. (C) Spearman's p correlation coefficient was calculated for the 
correlation between rank orders of every parameter. (D) Graphs represent the normalized ratio calculated as described in the materials and methods 
section. The black lines delimit the zone of compounds inducing a change inferior to DMSO ± 2aDMso> the red lines the zone of compounds 
inducing a change superior to DMSO ± 3aDMso- (E) Graphs represent the quantile-quantile (QQ) plots, plotting the theoretical quantiles for each assay 
in X-axis over the normalized ratio in y-axis. The grey points correspond to compounds inducing a change inferior to DMSO ± 2aDMso> the black points 
to a change comprised between DMSO ± 2aDMso DMSO ± 3aDMso> the red points to a change superior to DMSO ± 3aDMso. 



that the imidazole toxicity is independent of cholesterol depletion. 
Instead, imidazoles have an acute inhibitory effect on respiration, in 
line with previous reports^^. A potential explanation for this effect is 
the inhibition of the transient receptor potential cation channel, 
subfamily M, member 2 (TRPM2), a calcium-permeable cation 
channel that is regulated by free intracellular ADP-ribose^^. 
Importantly, even though the imidazoles are toxic in our experi- 
mental setup, they are mostly used topically in the clinic. In such 
application, less than 1% reaches the systemic circulation^^, which 
explains why this is unlikely to cause major systemic side effects. 

Impact of statins on metabolism in cultured cells and C. elegans. 

Similar to our strategy for imidazoles, we first performed dose- 
response studies for redox potential, ATP level and A\|/m using the 
different statins (Figure 6A-C). All three statins induced A\|/m by up 
to 3 -fold in Hepal.6 cells, indicating a strong hyperpolarization, 
while ATP level and redox potential did not change (Figure 6A- 
C). The effect on A\|/m was not due to an increase in the number 



of cells, which was only mildly increased at the highest dose with 
atorvastatin and simvastatin (Figure 6D). Treatment of Hepal.6 cells 
for 24 hours with fluvastatin and simvastatin decreased both OCR 
and ECAR, but only OCR with atorvastatin (Figure 6E-F), 
suggesting an overall decrease in metabolic rate of the cells. The 
decrease in OCR was maintained when OXPHOS was uncoupled 
(Figure 6G), indicating a decrease in overall respiratory capacity. 
However, unlike imidazoles, statins did not inhibit respiration in 
Hepal.6 cells acutely (Figure 6H). 

The effect of statins was further investigated in C. elegans using 
fluvastatin. This statin was selected as simvastatin and atorvastatin, 
but not fluvastatin, crystallized at 20 jiM in the agar plates used for 
worm culture. Although high doses of statins (mM) lead to embry- 
onic lethality in worms due to the impaired prenylation of small 
GTPases^^, we found that fluvastatin did not impair larval develop- 
ment up to the concentration of 50 |iM (Figure 61). We then tested 
the effect of fluvastatin on worm respiration (Figure 6J-K). While 
basal OCR was similar to untreated worms, addition of FCCP did not 
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Figure 4 | Hierarchical clustering and profiles of the 53 hits. (A) Hierarchical clustering of the 3 parameters for the 53 hits. Compounds were considered 
as hits when inducing a change superior to DMSO ± 2aDMso and when deviating from normal distribution (see methods section). Heatmap represents 
intensities of normalized ratios. (B-D) Remarkable profiles, i.e. combinations of the scores in each of the assays. (B) Profiles of toxic compounds that 
reduce significantly the three parameters. (C) Profiles of the cluster of compounds that increase significantly redox potential only. (D) Profiles of the 
cluster of compounds that increase significantly Av|/m only. 



increase OCR in fluvastatin -treated worms (Figure 6J-K), suggesting 
that the impact of statins on respiration was conserved through the 
different models. Finally, fluvastatin increased worm mean lifespan 
by 14% at 20 |iM and by 13% at 50 |iM (Figure 6L), indicating that 
the decrease in respiration induced by fluvastatin has an overall 
beneficial effect in worms. 

A potential mechanism to explain the decreased respiratory capa- 
city in cells and worms after long-term statin treatment is the inhibi- 
tion of mevalonate synthesis, which in turn limits CoQlO and 
heme-A levels^^, two essential components of the mitochondrial res- 
piratory chain^l Another possibility is that these effects are mediated 



at the transcriptional level. To test this hypothesis, we analyzed a 
dataset from the Connectivity Map database (http://lincs.hms. 
harvard.edu/db/). This data set contains gene expression profiles 
of the human breast cancer cell line MCF-7 exposed for 6 hours to 
three different statins, i.e. fluvastatin, simvastatin and lovastatin at 
10 |iM. The change in gene expression after statin treatment was 
low because of the limited number of samples per condition 
(Figure 7A). As expected, genes involved in cholesterol biosynthesis 
were strongly upregulated (Figure 7B). Additionally, we noted a 
significant downregulation of genes involved in mtDNA transcrip- 
tion (Figure 7C). 



SCIENTIFIC REPORT! | 4 : 5285 | DOI: 1 0.1 038/srep05285 



5 



ATP 



3 10 30 
concentration (|jM) 



DMSO 




1 3 10 30 
concentration (ijIVI) 

10|jM 



3 10 30 
concentration (iJIVI) 

20|jM 



3 10 30 
concentration (ijlVI) 

50mM 



s 













Basal ECAR 



Basal OCR 



3 10 30 
concentration (ijIVI) 



DMSO 



3 10 30 
concentration (ijIVI) 

10|jM 



G 

Uncoupled OCR 



H 



3 10 30 
concentration (ijlVI) 

20|jM 



Acute effect 
on OCR 



1.0- 
0.5- 
0.0- 
1.0- 
0.5 
0.0- 

1.0j».,._ 

0.5- 



0.0- 



10 100 
concentration (ijlVI) 

50|jM 



Figure 5 | Impact of imidazoles on cellular and C elegans metabolism. ( A-G) Effect of 24 hours treatment with econazole, miconazole, tioconazole and 
bifonazole in HepaL6 on ATP level (A), redox potential (B), Av|/m (C), DNA content (D), extra- cellular acidification rate (ECAR) (E), basal oxygen 
consumption rate (OCR) (F) and uncoupled OCR (G). (H) Acute effect of imidazoles on OCR in Hepal.6. Measurement was performed 3 minutes after 
addition of the compounds. Results are expressed as ratio over DMSO control and as mean ± SEM. Grey zone represents the DMSO 95% confidence 
interval. Significance was tested using a two-way ANOVA test, with Bonferroni correction for multiple comparison. * P<0.05, P<0.01 and 
P<0.001. Graphs represent mean ± SEM. (I) Econazole, miconazole, tioconazole and bifonazole impair the development of C. elegans N2 worms at all 
tested concentrations. Pictures were taken at day 3 of adulthood. 



As a further translation of our data towards human relevance, we 
performed a gene-set enrichment analysis on gene expression data 
from 480 patient- derived lymphoblastoid cell lines that were exposed 
for 24 hours to vehicle or 3 [iM simvastatin^^. The overall gene 
expression variation in these cells after simvastatin treatment was 
low, but due to the high number of samples per conditions, 26% of 
genes presented a significant change (Figure 7D). Similar to the 
Connectivity Map data, the cholesterol biosynthesis pathway was 
strongly induced by simvastatin (Figure 7E)^^'^^, and several mito- 
chondrial gene sets related to mtDNA transcription were downre- 
gulated (Figure 7F). In addition, expression of genes related to 
respiratory complex I and V were also repressed (Figure 7F). These 
data suggest that the mitochondrial functional changes we observed 
in mouse hepatocytes and worms are conserved in humans and are 
common to at least fluvastatin, simvastatin and lovastatin. 

Discussion 

We developed a set of straight-forward assays in two different model 
systems to rapidly identify whether a compound may have potential 
therapeutic benefits at the level of mitochondria. 

The first model is an immortalized hepatic cell line, Hepal.6. 
Immortalized cell lines offer the advantage that they can be grown 



in large quantities and high number of passages, two features that 
make their use easier when working at a high-throughput level. 
However, the immortalization process is accompanied by metabolic 
changes similar to those observed in cancer cell lines, including the 
presence of the so-called Warburg effect^^. A careful optimization of 
cell culture conditions is therefore required to ensure that the cells 
rely on mitochondrial metabolism for energy production. 
Supplementing oleic acid in the culture medium and increasing cell 
confluence stimulated Hepal.6 cells to switch from glycolytic to 
mitochondrial metabolism (Figure 2). In this cell line, we measured 
several parameters that displayed a wide dynamic range during our 
screening. This allowed the identification of both activators and 
inhibitors of mitochondrial function. In fact, we were able to stratify 
the hits according to their response and divided into them separate 
classes, i.e. compounds toxic to mitochondria and cells (e.g. papa- 
verine, niclosamide), compounds that induce redox potential (imi- 
dazoles) and those that induce A\|/m signal (statins). 

Hits from our high-throughput screening were selected for the 
evaluation of respiration and lactate production. The increase in 
the throughput of the respiration assays enabled the measurement 
of up to four parameters related to catabolism (ECAR, basal OCR 
and uncoupled OCR after long-term treatment, and OCR after acute 
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Figure 6 | Impact of statins on cellular and C elegans metabolism. (A-G) Effect of 24 hours treatment with atorvastatin, fluvastatin and simvastatin in 
Hepal.6 on ATP level (A), redox potential (B), Av|/m (C), DNA content (D), extra- cellular acidification rate (ECAR) (E), basal oxygen consumption rate 
(OCR) (F) and uncoupled OCR (G). (H) Acute effect of statins on OCR in Hepal.6. Measurement was performed 3 minutes after addition of the 
compounds. (I) Fluvastatin does not impair the development of C. elegans^l worms at 10, 20 and 50 [iM. Pictures were taken at day 3 of adulthood. (J- 
K) Effect of fluvastatin on basal OCR (J) and uncoupled OCR (K) at day 1 of adulthood. (L) Lifespan in N2 worms after treatment with fluvastatin at 20 
and 50 [iM. Results are expressed as ratio over DMSO control and as mean ± SEM. Grey zone represents the DMSO 95% confidence interval. Significance 
was tested using a two-way ANOVA test, with Bonferroni correction for multiple comparison. * P<0.05, P<0.01 and *** P<0.001. 



treatment). Notably, the determination of OCR after acute incuba- 
tion with the compounds allowed us to distinguish between the direct 
impact of imidazoles on mitochondria versus the indirect effect of the 
statins on respiration, only evident after long-term treatment. 

The second model that we used to analyze mitochondrial function 
was the roundworm C. elegans. It has the advantage to provide 
insights on the effects of compounds in an intact animal and at the 
scale of multiple tissues. One of our read- outs was the impact of 
compounds on larval development. C. elegans development depends 
on the accurate timing of cell division and differentiation, which 
requires a complex combination of genetic and molecular mechan- 
isms^^. In particular, development of C.elegans is highly energy- 
dependent and relies on the electron transport chain (ETC), as it 
was shown that food deprivation and downregulation of ETC com- 
plexes expression strongly delayed or stopped the larval develop - 
ment^°. In line with these premises, imidazoles, but not statins, 
caused a severe developmental delay. Since C. elegans are cholesterol 
auxotrophs, our data suggest that the toxicity of the imidazoles is 
independent of cholesterol depletion. Instead, imidazoles have an 
acute inhibitory effect on respiration, which maybe due to the induc- 
tion of uncoupling and swefling of mitochondria, as previously docu- 
mented in isolated rat liver mitochondria exposed to econazole^^. 

In contrast, the mechanism by which statins affect mitochondrial 
respiration was not associated with toxicity, as Hepal.6 cefls were 
able to maintain their ATP and redox potential levels, and worm 
development was not impacted. Our result hence contrasts with 
previous reports showing that high pharmacological doses of statins 
inhibit C. elegans development^'*. We used, however, at least ten times 
lower doses, and stayed below doses at which the hydrophobic flu- 
vastatin crystaflizes in agar plates. Working in these more gentle and 
physiological conditions, we were able to show for the first time that 
fluvastatin increased worm lifespan, in line with the lifespan exten- 



sion observed in Drosophila and the reduction in afl-cause mortality 
measured in humans after simvastatin treatment^*"^^. Remarkably, 
the fact that statins extend lifespan in lower organisms suggests that 
their therapeutic benefit in human may also be partly independent of 
their effect on cholesterol biosynthesis. 

Two potential mechanisms can explain the decreased respiratory 
capacity in cefls and worms after long-term statin treatment. First, by 
inhibiting mevalonate synthesis, statins also impact on the levels of 
CoQlO and heme-A^^, which are both essential for respiratory chain 
function^^ Second, gene expression data from both lymphoblastoid 
cefl lines and MCF-7 cefls show that there is also a downregulation of 
genes involved in mtDNA transcription. In the lymphoblastoid cell 
lines, there is furthermore a pronounced decrease in the expression 
of genes encoding for components of complex I and V subunits. The 
link between these data and the increase in A\|/m after statins is less 
clear. In fact, A\\f^ is decreased in cefls with defects in respiratory 
complexes I and IP^'^^. However, another study showed that in neu- 
rons deficient for respiratory complex I, Ay\f^ was increased because 
ATP synthase is working in reverse mode^^. Overafl, these data 
underline the complexity of the metabolic remodeling induced by 
statins treatment. This complexity is directly linked to the number of 
pathways affected by the inhibition of mevalonate synthesis, which is 
an early precursor of multiple metabolites. 

In conclusion, we developed a cross-species mitochondrial screen- 
ing and validation platform, using both mammalian cell lines and C. 
elegans, exploiting and emphasizing the strengths of each of these 
models. Using this screen, we identified negative and positive regu- 
lators of mitochondrial function. The combination of multiple tests 
that inform about distinct aspects of mitochondria enabled us to 
define a complex and comprehensive signature of the actions of 
pharmacological compounds on mitochondria. This setup allowed 
us to define the effects of imidazoles and statins on mitochondrial 
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Figure 7 | Impact of simvastatin on gene expression in MCF-7 cells and human lymphoblastoid cell lines. (A-C) Effect of three statins, fluvastatin (Flu.), 
lovastatin (Lov.) and simvastatin (Sim.), on gene expression in MCF-7 cells after 6 hours treatment at 10 \iM. (A) Vulcano plot of microarray dataset. 
Treatment/control fold change expression was determined by calculating first the average expression for all statins and DMSO treatments. Cholesterol 
synthesis geneset (GO:0006695) and mitochondria geneset (GO:0005739) are represented in red and blue, respectively. Black horizontal bar corresponds 
to the significance threshold for a nominal p-value of 0.05, while the red horizontal bar corresponds to the p-value for multiple testing after Bonferroni 
correction. (B) Positive enrichment score indicates that MCF-7 cells induce cholesterol synthesis genes (geneset GO:0006695, p<0.01) following sterol 
depletion. (C) Negative enrichment score for mtDNA transcription genes (geneset GO:0006390, q=0.24) shows that mtDNA processing is 
downregulated after statins treatment in MCF-7 cells. (D-F) Effect of simvastatin on gene expression in lymphoblastoid cell lines derived from patients 
after 24 hours treatment at 3 [iM. (D) Vulcano plot of microarray dataset. Cholesterol synthesis geneset (GO:0006695) and mitochondria geneset 
(GO:0005739) are represented in red and blue, respectively. The red horizontal bar corresponds to the p-value for multiple testing after Bonferroni 
correction. (E) Positive enrichment score indicates that in lymphoblastoid cells simvastatin induces cholesterol synthesis genes (geneset GO:0006695, 
p<0.001). (F) Negative enrichment scores for mtDNA transcription genes (genesets GO:0006264, q=0.009; GO:0032042, q=0.003 and GO:1901858, 
q=0.084), respiratory complex I (genesets GO:0005747, q=0.052 and GO:0032981, q=0.143) and respiratory complex V (genesets GO:005753, q=0.015; 
GO:0042776, q= 0.039 and GO:0070071, q= 0.079) show that respiratory complexes and mtDNA processing is downregulated after simvastatin treatment 
in lymphoblastoid cells. For B, C, E and F, genesets were determined according to Gene Ontology annotation (www.geneontology.org). Heatmaps 
represented in E and F show only the genes that are significantly differentially expressed between control (green bar) and statin (red bar) treated samples 
according to nominal p values. 



function, and identify potential novel pathways that may contribute 
to the pleiotropic effects observed upon statin treatment. We envi- 
sion that such a comprehensive screening approach is essential to 
identify novel mitochondrial modulators. Such compounds will be 
beneficial for both rare and common diseases in which mitochondria 
are dysfunctional. 

Methods 

Compounds used in cell-based assays and C. elegans. The Prestwick Chemical 
Library, which contains 1,200 small molecules and FDA approved drugs, was 
acquired at Prestwick Chemicals (lllkirch-Graffenstaden, France). Fluvastatin, 
simvastatin and atorvastatin were provided by Prestwick Chemicals (lllkirch- 
Graffenstaden, France). Econazole, miconazole, tioconazole and bifonazole were 
acquired at Sigma (Sigma- Aldrich, St. Louis, MO, USA). 



Cell lines and medium conditions. Hepal.6 mouse hepatoma cells were obtained 
from ATCC and maintained according to ATCC guidelines. For the assays, hepal.6 
cells were plated in low glucose medium (1 g/1 glucose) supplemented with 30 [lM 
oleate-BSA (LGO) (Sigma- Aldrich, St. Louis, MO, USA), unless specified otherwise. 

High-throughput cell-based assays. Hepal.6 cells were seeded at 8,000 cells per well 
in a volume of 15 |il, in 384-well plates. Compounds and were added from a 4x 
concentrated stock, diluted in culture medium, to the cells, reaching a total assay 
volume of 20 |il. Assays were typically performed 24 hours after incubation with 
medium or compounds. ATP levels were measured using Cell Titer Glo (Promega, 
Madison, WI, USA). Mitochondrial membrane potential (A\|/m) was measured using 
the fluorescent tetramethylrhodamine methyl ester (TMRM) probe (T-668, 
Invitrogen, Paisley, UK). Briefly, cells were incubated in culture medium and TMRM 
was added to the culture medium at a final concentration of 200 nM for 45 minutes at 
37°C before washing with 50 [l\ of PBS. The Alamar Blue assay (Life Technologies, 
Paisley, UK) was used as a marker for redox potential, following manufacturer 
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protocols. ATP level, redox potential and TMRM assays were measured on a Tecan 
Infinite 500 (Tecan, Mannedorf, Switzerland), and performed in duplicate. 

Screening of a 1,200 drugs library - specific procedures. Hepal.6 cells were seeded 
at 8,000 cells per well, in a volume of 20 in 384-well plates, treated on day 1 using a 
Sciclone G3 Liquid Handling Workstation (PerkinElmer, Waltham, MA, USA) and 
distributing a volume of 2.22 |il of compounds at 100 |iM to obtain a final 
concentration of 10 |iM in each well. Assays were performed after 24 hours of 
incubation with compounds. For ATP level, redox potential and TMRM assays, we 
used DMSO at 0.1% as a negative control (final DMSO concentration in all wells) and 
oligomycin A at 1 |J,M as a positive control, both of them being present on each side of 
all the plates, i.e. columns 23-24 and 1-2 respectively. These assays were performed in 
duplicate. 

High-throughput screening data were analyzed using the cellHTS2 package from R 
software (http://www.r-project.org/). ATP, redox potential and TMRM, data were 
normalized per plate using the normalized percentage of inhibition function (nor- 

pos 

malized ratio), calculated as follows: x^f^ = — ^ , where x^i is a result point, /if 

is the average for all the oligomycin A positive control points and fi^^^ is the average of 
the DMSO negative control points. Z scores and Z factors were calculated from these 
normalized values for the Figure 3 A. Theoretical quantiles (Figure 3E) were calculated 
from the normalized values using R software. Compounds were considered as hits 
when inducing a signal change superior to DMSO ± 2g^mso and when deviating 
from normal distribution, i.e. when squared distance between observed quantile and 
theoretical quantile was above a certain threshold. This threshold was arbitrarily set at 
> 0.005 for ATP, TMRM and redox potential signals. Clustering of the 53 hits was 
performed on the Hierarchical Clustering module available on the GenePattern 
website (http://www.broadinstitute.org/cancer/software/genepattern), using 
Euclidean distance and single linkage as analysis parameters. 

Measurement of cellular NAD^ and NADH levels. NAD^ and NADH levels were 
determined using the EnzyChrom™ NAD/NADH Assay kit and following 
manufacturer's instructions (BioAssays Systems, Hayward, CA, USA) 

Measurement of DNA levels per cell. DNA was measured directly in Seahorse 96- 
well plates using the CyQUANT® Cell Proliferation Assay Kit and following 
manufacturer's intructions (Lifetechnologies, ThermoFisher Scientific, Waltham, 
MA, USA). 

Measurement of cellular and C. elegans respiration. Oxygen consumption rate 
(OCR) and extra- cellular acidification rate (ECAR) were measured using the 
Seahorse XF96 equipment (Seahorse bioscience Inc., North Billerica, MA, USA). For 
measurement of cellular metabolism, Hepal.6 were seeded at 12,000 cells per well and 
treated with different medium conditions or compounds for 24 hours in 150 |il 
medium. Acute effect of compounds was evaluated by injecting them during the run. 
Worm oxygen consumption was performed as described''*^. Briefly, at day 1 of 
adulthood, ten worms per well were transferred into M9-filled 96-well Seahorse 
plates. Oxygen consumption was measured 6 times in basal conditions. FCCP 
treatment was performed at 10 [lM final. 

Pharmacological treatment of C. elegans strains. C. elegans strains were cultured at 
20°C on nematode growth media agar plates seeded with E. coli strain HTl 15 (DE3). 
Strains used were wild-type Bristol N2 and provided by the Caenorhabditis Genetics 
Center (University of Minnesota). Imidazoles and statins were added at the indicated 
concentrations just before pouring the plates. Animals were exposed to compounds 
during the full life from hatching onwards. To ensure a permanent exposure to the 
compound, plates were changed twice a week for the lifespan assays. Imidazoles and 
statins were dissolved in DMSO, and the vehicle controls were adapted accordingly. 

Lifespan and development assays. Lifespan tests were performed as described at 
20°G*^. Briefly, 60-100 animals were used per conditions and scored every other day. 
Animals that crawled off the plate or had an «exploded vulva» phenotype were 
censored. The development of a synchronized worm population was monitored after 
4 days of post-embryonic development by taking pictures of the whole plate. 

Worm statistics. Survival analyses were performed using the Kaplan Meier method 
and the significance of differences between survival curves calculated using the log 
rank test. Calculation of mean lifespan and SEM were calculated using the R software. 

Bioinformatics analysis. The datasets generated in human lymphoblastoid cell 
lines^^ was downloaded from Gene Expression Omnibus under the accession number 
GSE36868. The datasets generated in the human breast cancer cell line MCF-7 was 
downloaded from the Connectivity Map database (http://lincs.hms.harvard.edu/db/). 
Data were first normalized with RMAExpress software (http://rmaexpress. 
bmbolstad.com/) and analyzed using the Gene Set Enrichment Analysis software^". 
The cholesterol synthesis geneset was defined according to GO annotation 
(Accession number: 0006695). In this case, since a single gene set was tested (e.g. 
cholesterol gene set), a nominal p value inferior to 0.05 was used as significance 
threshold. Mitochondrial genesets were defined by searching for all gene ontologies 
containing the word "mitochondria*" on AmiGO website (http://amigo. 
geneontology.org/cgi-bin/amigo/go.cgi). A total of 121 genesets were defined using 



this method. In this case, several genesets were tested and the q value of the false 
discovery rate control inferior to 0.25 was used as a significance threshold. 
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